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In  the  near-infrared  and  visible  bandpasses  optical  propagation  theory  conventionally  assumes  that  hu¬ 
midity  does  not  contribute  to  the  effects  of  atmospheric  turbulence  on  optical  beams.  While  this  assumption 
may  be  reasonable  for  dry  locations,  we  demonstrate  that  there  is  an  unequivocal  effect  owing  to  the 
presence  of  humidity  upon  the  strength  of  turbulence  parameter,  C„2,  from  data  collected  in  the  Chesa¬ 
peake  Bay  area  over  100  m  length  horizontal  propagation  paths.  We  describe  and  apply  a  novel  technique, 
Hilbert  phase  analysis,  to  the  relative  humidity,  temperature,  and  C„2  data  to  show  the  contribution  of 
the  relevant  climate  variable  to  C„2  as  a  function  of  time.  ©  2007  Optical  Society  of  America 
OCIS  codes:  010.1300,  010.7060,  030.7060. 


1.  Introduction 

It  has  been  known  for  some  time1  that  the  scintilla¬ 
tion  behavior  of  point  sources  is  a  measure  of  the 
optical  seeing  in  the  atmosphere.  What  has  been  less 
well  understood  is  the  contribution  of  different  envi¬ 
ronmental  variables  to  optical  seeing. 

Over  the  past  decade,  a  great  deal  of  study  has  been 
dedicated  to  clarifying  this  issue.  Comprehensive 
treatments  of  the  theory  of  wave  propagation  in  ran¬ 
dom  media  are  given  in  Tatarskii’s  seminal  works.2’3 
Some  of  the  simplest  models  based  on  these  complex 
works  are  well  known  and  available  in  the  literature: 
Hufnagel-Valley,4  SLC-Day  and  SLC-Night.5  These 
models  are  commonly  used  to  predict  the  strength  of 
weak  clear  air  turbulence’s  refractive  index  structure 
parameter,  C2. 

The  underlying  assumption  is  that  the  index  of  re¬ 
fraction  of  air  depends  solely  on  pressure  and  temper¬ 
ature  at  visible  and  near-infrared  wavelengths.2’6  We 
can  verify  this  from  landbound  horizontal  path  mea¬ 
surements.  These  show  that,  during  daylight  hours  on 
a  clear  sunny  day,  the  weak  scintillation  regime  C  2 
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trend  is  dominated  by  solar  insolation,  as  illustrated  in 
the  left-hand  plot  of  Fig.  1. 

The  right-hand  plot  of  Fig.  1  illustrates  the  well- 
known  fact  that  the  effects  of  turbulence  do  not  cease 
after  sunset.  Here  it  is  less  clear  as  to  the  predominant 
contributing  factors.  The  assumption  of  temperature 
dominance  is  sensible  in  the  driest  of  environments, 
but  we  will  demonstrate  from  experimental  evidence 
that  this  is  an  incomplete  condition  where  there  is  a 
significant  amount  of  water  vapor  entering  or  leaving 
air  volume  where  G2  is  measured. 

A.  Previous  Work 

A  body  of  work  has  been  reported  whose  aims  have 
been  to  compare  bulk  climate  estimates  of  C2  with 
optically  derived  readings.7-9  In  those  works,  the 
wavelengths  of  the  optical  measurements  were  in  the 
mid-infrared  (3  to  5  pm),  and  all  showed  that  humid¬ 
ity  along  the  measurement  volume  plays  a  significant 
role  in  the  temporal  behavior  of  C2. 

In  a  theoretical  study  of  the  structure  parameter 
using  bulk  measurables  over  sea  ice  and  snow,  An¬ 
dreas10  defined  a  refractive  index  scale  n*  via 

n*  Bq* 

n*=At*  +  Bq*  -a  ^  =  1  +  (1) 

where  t,  and  q,  represent  temperature  and  humidity, 
scaled  such  they  are  constant  with  height  and  A  and 
B  are  constants  given  in  Appendix  A.  The  variable  n  , 
therefore  represents  a  horizontal  refractive  index 
layer  above  the  terrestrial  surface.  He  rewrote  Eq.  (1) 
in  terms  of  the  Bowen  ratio  Bo,  the  ratio  of  sensible 
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Fig.  1.  Effects  of  solar  insolation  on  C„2  as  measured  on  9  March 
2006  in  Puerto  Rico.  Solar  radiation  measurements  are  superim¬ 
posed  at  an  arbitrary  scale  on  top  of  the  C„2  data.  (Left)  The  full 
24  h  period.  The  vertical  axis  values  are  of  order  10~12nT2/3.  (Right) 
Close  up  of  the  evening  period.  The  vertical  axis  values  are  of  order 
10~14irT2/3.  The  C„2  data  were  obtained  at  an  urban  site  1.75  km 
from  the  sea.  These  values  are  in  agreement  with  measurements  of 
seeing  over  the  sea  by  alternative  means  (Ref.  23). 

heat  flux  to  latent  heat  flux,  such  that 

n*  /p  cp\  B  j  B  \ 

AF„  =  1  +  ( X  )A(Bo )  =  1  +  \KA(Bo ) )’  (2) 

where  cp  represents  the  specific  heat  of  air  at  constant 
pressure,  L  is  the  latent  heat  of  the  sublimation  of  ice, 
and  p  is  the  air  density.  Since  Bo  °c  t,/qt,  this  formu¬ 
lation  allows  the  dependence  of  the  refractive  index 
scale  on  temperature  and  humidity  scales  to  be  easily 
interpreted. 

We  adapt  Eq.  (2)  to  our  experimental  conditions  by 
setting  L  to  be  the  latent  heat  of  vaporization  of  water, 
and  we  assume  that  P  =  1000  hPa,  T  =  25  °C,  and 
X  =  0.93  pm.  The  result  is  shown  in  Fig.  2;  from  this 
it  is  clear  that,  when  \Bo  |  is  large,  n,/At ,  is  near  1. 
This  means  that,  when  sensible  heat  (heat  exchange 
without  change  of  thermodynamic  phase)  dominates 
the  Bowen  ratio,  the  t,  in  the  denominator  of  Eq.  ( 1) 
is  the  principal  parameter. 

When  |  Bo  |  is  small,  however,  the  latent  heat  flux 
dominates,  and  q,  makes  the  major  contribution  to  n,. 
Thus  there  is  an  interplay  between  humidity  and 


temperature  effects  on  the  refractive  index  in  the 
visible  and  near  infrared,  which  varies  as  a  function 
of  the  Bowen  ratio.  The  pole  at  the  origin  of  Fig.  2  is 
indicative  of  a  breakdown  of  this  formulation  at  very 
small  Bowen  ratios.  Nevertheless  the  curve  strongly 
suggests  that,  as  long  as  there  is  significant  injection 
of  moisture  into  the  measurement  volume  of  C2,  the 
humidity  effects  become  nonnegligible. 

Here  we  demonstrate  that  the  temporal  behavior  of 
the  optically  determined  turbulence  structure  param¬ 
eter  C  2  in  the  visible  and  near-infrared  bandpasses  at 
a  coastal  site  is  strongly  influenced  by  local  humidity. 
Although  no  Bowen  ratio  values  could  be  determined 
from  the  available  measurements,  the  results  do  sup¬ 
port  the  qualitative  interpretation  of  Andreas’s  work. 

2.  Experiment  and  Correlogram  Analysis 

The  C2  and  associated  weather  variable  data  used  in 
this  study  were  collected  over  a  number  of  days  during 
2003  and  2004  at  Chesapeake  Beach  in  Maryland, 
USA.  The  C2  data  were  obtained  with  a  commercially 
available  scintillometer  [Model  LOA-004,  Optical  Sci¬ 
entific,  Gaithersburg,  Maryland] ,  which  served  as  both 
a  scintillometer  and  as  an  optical  anemometer  for 
winds  transverse  to  the  beam  paths.  The  separation 
between  source  and  receiver  was  just  under  100  m, 
placing  the  system  in  the  weak  scintillation  regime. 
The  local  weather  parameters  were  determined  by  a 
Davis  Provantage  Plus  (DP+)  weather  station.  The 
LOA-004  had  a  sample  rate  of  10  s,  and  the  DP+ 
was  set  at  5  min. 

The  LOA-004  comprises  a  single  modulated  in¬ 
frared  transmitter  whose  output  is  detected  by  two 
single  pixel  detectors.  The  transmitting  LED  has 
a  bandwidth  of  0.92  to  0.96  pm,  while  the  detec¬ 
tor  bandwidth  is  much  broader,  at  0.65  to  1.0  pm. 
The  path-integrated  C2  measurements  are  deter¬ 
mined  by  the  LOA-004  by  computation  from  the 
log-amplitude  scintillation  [Cx(r)j  of  the  two  receiv¬ 
ing  signals.11’12  The  algorithm  for  relating  Cx(r)  to  C2 
is  based  on  an  equation  for  the  log-amplitude  covari¬ 
ance  function  in  Kolmogorov  turbulence  by  Clifford 
et  al.13 

The  C2  data  were  smoothed  with  a  60  point  (5  min) 
rolling  average  function.  We  define  the  morning  and 
night  portions  of  a  24  h  period  with  respect  to  the 
measured  solar  irradiance  function  such  that  we  ex¬ 
clude  the  effect  of  solar  insolation  from  the  data  in 
this  study.  Morning  runs  from  midnight  until  sun¬ 
rise,  while  night  runs  from  sunset  until  23:59.  As 
reported  in  Oh  et  al.14~16  visual  inspection  of  the  time 
series  data  gives  the  impression  that  there  is  an  ap¬ 
proximate  inverse  relationship  between  C  2  and  rel¬ 
ative  humidity.  This  can  be  appreciated  in  a  more 
quantitative  manner  by  graphing  one  parameter 
against  the  other. 

We  chose  data  sets  in  which  the  temperature  vari¬ 
ations  are  no  more  than  16  °F,  and  the  pressure 
change  is  at  most  15  hPa  over  the  time  intervals  of 
interest.  The  actual  range  of  variation  per  data  set  is 
shown  in  Table  1.  The  data  sections  were  also  se- 
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Table  1.  Mean  and  Range  of  Bulk  Parameter  Measurements  and 
the  Maximum  Range  of  Specific  Humidity  (A0s) 


Date 

RH 

(%) 

T 

(°F) 

P 

(hPa) 

A RH 

\T 

Ap 

A Qs 

(g/kg) 

11/03/03  (p.m.) 

83.6 

65.2 

1015.3 

15 

2 

0.4 

2.2 

11/09/03  (a.m.) 

49.3 

38.8 

1030.9 

11 

3.1 

1.2 

0.5 

11/10/03  (a.m.) 

76.4 

46.8 

1027.6 

11 

8 

1.8 

0.6 

02/02/04  (a.m.) 

91.4 

29.9 

1022.9 

16 

2.5 

6.6 

0.8 

03/27/04  (p.m.) 

76.4 

46.8 

1027.6 

11 

8 

1.8 

1.7 

03/28/04  (a.m.) 

71 

44.5 

1027.6 

16 

0.7 

0.6 

0.9 

04/03/04  (a.m.) 

76.1 

46.6 

997.4 

18 

5.9 

4.9 

0.4 

04/03/04  (p.m.) 

43 

40.1 

1004.5 

14 

5.6 

2.1 

1.1 

lected  to  have  no  scattering  effects  attributable  to 
snow  or  rain,  and  the  wind  was  northerly  (to  within 
approximately  ±20°,  inflowing  from  the  bay  to  land). 
Only  a  small  subset  of  eight  morning  and  evening 
runs,  spanning  seven  days  between  November  2003 
and  March  2004,  provided  a  complete  time  series  in 
both  ambient  weather  variables  and  C,,2.  Part  of  this 
is  shown  in  Fig.  3;  the  Cn2  versus  humidity  correlo- 
grams  evidence  a  negative  gradient  in  all  eight  runs. 
Table  1  shows  us  that  the  relative  humidity  variation 
is  more  strongly  reflective  of  changes  in  absolute 
moisture  content  than  of  moisture  holding  capacity  of 
the  air. 

We  conclude  that  humidity  plays  a  significant  part 
in  the  behavior  of  the  refractive  index  structure  pa¬ 
rameter  C„2  in  the  data  sets  studied. 


3.  Stationarity  Problem 

The  unsophisticated  use  of  the  correlogram  tool  in 
Section  2  is  a  rapid,  first-order  method  for  inferring 
the  statistical  influence  of  one  measurable  upon  an¬ 
other.  However  it  is  unsatisfactory  as  it  does  not 
reveal  any  detailed  information  such  as  exactly  when 
the  humidity  contribution  is  important  with  respect 
to  all  other  possible  parameters  (e.g.,  temperature) 
and  to  what  degree  its  contribution  is  influential. 
Cross-covariance  techniques  are  a  possible  second- 
order  method17  for  extending  the  analysis,  but  there 
are  two  major  stumbling  blocks.  The  first  is  the  dif¬ 
ficulty  in  interpreting  the  meaning  of  the  covariance 
functions  physically,  and  the  second  is  the  nonsta¬ 
tionary  nature  of  the  data.  In  the  latter  case  it  is  well 
known  that  strict  stationarity  is  a  constraint  that  is 
impossible  to  satisfy  on  practical  grounds  since  no 
detector  can  cover  all  possible  points  in  phase  space. 
This  has  motivated  us  to  employ  a  novel  analysis 
technique  to  infer  the  relationship  between  C,,2  and 
the  climate  parameters,  which  we  briefly  describe 
in  the  following  section.  We  leave  a  more  detailed 
study  of  the  limitations  of  the  technique  to  a  sepa¬ 
rate  paper. 

A.  Analytic  Signal 

Gabor18  defined  the  complex  analytic  signal,  namely, 
^(t)  =X(t)  +  iY(t), 

where 


03/27/2004  p.m.  (T  =46.8"  F)  03/28/2004  a.m.  (T  =44.5°  F) 

r  '  av  av 


04/03/2004  a.m.  (T  =46.6  F)  04/03/2004  p.m.  (T  =40.1°  F) 

'  av  '  r  '  av  ' 


Fig.  3.  Example  correlograms  of  C„2  and  relative  humidity  in  the 
absence  of  solar  insolation.  The  upper  and  lower  bounds  indicate 
the  50%  confidence  level.  The  C„2  magnitudes  are  10~15m~2/3,  in 
agreement  with  measurements  of  seeing  over  the  sea  by  alterna¬ 
tive  means  (Ref.  23). 


1  f  X(-) 

Y(t)=nx(t)]=  ^  pj  (^)dT’  <3> 


where  W[*]  represents  a  Hilbert  transform.  Through 
this  definition,  'V(t)  is  unique.  The  Hilbert  transform 
is  a  well-known  integral  transform  with  a  singular 
kernel  (1/[tt (t  -  t)]),  t  also  being  a  time  variable  if  t 
is  time.  It  is  also  a  Cauchy  principal  integral,  which 
we  denote  by  P.  Note  that  the  Hilbert  transform  pre¬ 
serves  the  norm  of  the  real  signal;  the  difference  be¬ 
tween  Y  and  X  is  a  tt/2  phase  shift  for  a  periodic 
function.  That  is,  the  Hilbert  transform  of  cos  t  is 
sin  t,  and  sin  t  is  —cos  t,  and  the  transform  of  a 
constant  is  zero. 

Equation  (3)  means  that  we  can  write  the  analytic 
signal  as 


M'(7)  =  a(t)  exp1 


M>(0 


where  a(t )  =  \X2(t)  +  Y2(t), 


=  arctan 


(m\ 

\mr 


(4) 


which  is  similar  to  the  well-known  Fourier  expres¬ 
sion.  We  may  now  determine  the  instantaneous  sig¬ 
nal  phase,  <!>(/).  Thus  we  can  also  calculate  the 
instantaneous  frequency  co(t),  defined  as 
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“(0 


d<b(f) 

dt 


(5) 


Both  d>(t)  and  w(f)  can  be  interpreted  as  physical 
measurables  provided  certain  preconditions,  which 
we  describe  below,  are  met. 

B.  Hilbert  Phase  Analysis 

There  is  a  problem  with  the  definition  of  Eq.  (3);  the 
Hilbert  transform’s  kernel  represents  a  noncausal  fil¬ 
ter,  of  infinite  support.  So  if  one  applies  the  Hilbert 
transform  directly  to  a  time-varying  signal  that  has 
a  nonzero  local  mean  in  any  subsection,  there  is  a 
high  probability  that  at  least  one  of  a  number  of 
paradoxes19  will  be  encountered.  This  may  be  best 
appreciated  if  we  consider  a  signal  in  phase  space,  as 
in  the  left-hand  plot  of  Fig.  4.  Here  we  show  the  real 
versus  the  imaginary  components  of  the  Hilbert 
transform  of  a  series  of  analytic  C„2  measurements. 
The  trajectory  of  the  analytic  signal’s  vector  is  subject 
to  many  alterations  in  the  vector  start  point,  norm, 
and  phase  angle.  Any  attempt  to  determine  the  in¬ 
stantaneous  frequency  is  bound  to  be  problematic  as 
we  are  not  able  to  follow  the  signal  vector’s  start 
position  over  time;  this  starting  point  obviously  does 
not  remain  at  the  origin  of  the  coordinate  system. 
Calculating  the  instantaneous  frequency  from  the 
Hilbert  transform  as  is  will  generate  both  positive 
and  negative  values,  rendering  it  physically  uninter¬ 
pretable. 

The  paradoxes  may  be  avoided  by  application  of 
the  empirical  mode  decomposition  (EMD)  method 
developed  by  Huang  et  al.,20  which  we  have  imple¬ 
mented.21-22  EMD  is  a  unique  and  novel  method 
that  is  able  to  separate  an  arbitrary  real-time  series 
into  eigenfunctions  termed  intrinsic  mode  functions 
or  IMFs,  each  of  which  possesses  a  structure  with 
well  defined  instantaneous  frequencies,  u)(t).  The 
term  eigenfunction  is  used  suggestively  here;  we  do 
not  mean  to  imply  that  the  IMFs  are  eigenfunctions 
in  the  strict  sense. 


Briefly,  the  EMD  technique  consists  of  (1)  the  de¬ 
termination  of  two  envelope  functions  about  the  time 
series  covering  all  the  local  minima  and  maxima, 
respectively;  (2)  the  computation  of  the  mean  of  the 
two  envelope  functions;  (3)  the  removal  of  the  mean  of 
the  envelopes  from  the  time  series  in  an  iterative 
manner  until  the  mean  is  found  to  be  zero;  and  (4) 
storing  of  the  resulting  mode  as  an  IMF.  The  IMF  is 
then  subtracted  from  the  original  time  series  and 
steps  (1)  to  (4)  are  repeated,  thereby  sifting  out  a 
family  of  IMFs,  stopping  only  when  the  resulting 
mode  shows  no  variation.  The  final  mode  represents 
the  overall  trend  of  the  signal  and  is  not  itself  an  IMF. 

The  Hilbert  transform  of  the  IMFs,  one  of  which  is 
shown  on  the  right-hand  side  of  Fig.  4,  ensures  that 
the  analytic  signal  vector’s  origin  stays  fixed  and  no 
sudden  changes  in  the  direction  of  oAt)  occur.  These 
conditions  being  satisfied,  the  instantaneous  fre¬ 
quency  remains  positive  and  physically  meaningful. 

4.  Hilbert  Phase  Analysis 

The  Hilbert  phase  analysis  (HPA)  technique  is  based 
on  the  ideas  mentioned  in  the  previous  subsections.  It 
is  clear  that  a  phase  angle  [<!>(£)],  as  well  as  an  ampli¬ 
tude  [a(D],  can  be  found  from  the  Hilbert  transform  of 
the  IMFs  derived  from  the  EMD  sifting  process. 

The  procedure  for  the  HPA  data  analysis  follows  a 
three-step  protocol: 

1.  Decompose  the  time  series  measurements  of 
differing  parameters  via  EMD,  obtaining  their  IMFs 
and  trend  lines. 

2.  Apply  the  Hilbert  transform  to  the  various 
IMF  sets.  We  discard  the  trend  lines  as  the  reason 
that  the  original  time  series  fail  to  be  Hilbert  trans¬ 
formed  in  a  physically  comprehensible  manner. 

3.  Examine  the  instantaneous  phase  angles  of 
the  Hilbert  transformed  IMFs  between  different 
parameters  to  infer  the  dynamics  of  the  physical 
system. 


Fig.  4.  (Left)  Hilbert  phase  space  plot  of  the  trajectory  of  the  C„2 
signal  vector  of  Fig.  1.  The  signal  vector’s  start  point  drifts  arbi¬ 
trarily  around  the  phase  space  making  the  determination  of  a 
physically  reasonable  instantaneous  frequency  impossible  over  the 
whole  path.  (Right)  Hilbert  phase  space  plot  of  the  trajectory  of  a 
single  IMF  derived  from  the  C„2  signal.  The  IMF  vector’s  start 
point  is  stable,  and  its  trajectory  does  not  change  direction,  so  a 
positive  instantaneous  frequency  can  be  determined  at  all  points. 


A.  Interpretation  of  Hilbert  Phase  Analysis 

Here  we  demonstrate  that  the  physical  effects  of  a 
nonlinear,  nonstationary,  time-varying  system  can 
be  studied  via  the  sum  of  all  IMF  phases,  2imf  dKt). 

In  Fig.  5  the  Simf  Th)  is  graphed  with  the  solar 
radiation  data  superimposed  for  example  days.  There 
is  clearly  a  gradient  change  in  the  phase  function  at 
sunrise  and  sunset.  Also  notable  is  the  phase  jump 
toward  a  lower  gradient  whenever  the  solar  radiation 
function  exhibits  a  drop  in  amplitude;  likewise  the 
phase  gradient  increases  with  sudden  increases  in 
the  measured  solar  radiation.  This  can  be  understood 
as  a  reduction  in  energy  in  the  system  leading  to  a 
lower  instantaneous  frequency,  therefore  we  see  a 
lower  phase  gradient.  A  change  in  energy  results  in 
a  change  in  instantaneous  frequency,  so  we  see  a 
modified  gradient. 

We  conclude  that  the  overall  physical  features  of  a 
nonstationary  time  series  can  be  extracted  by  an  in¬ 
spection  of  the  sum  of  its  Hilbert  phases  of  the  IMFs. 
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Fig.  5.  Sum  of  all  Hilbert  phases  of  the  measured  solar  radiation 
IMFs  f°r  the  following  days:  2  February,  27  March,  28 

March,  and  3  April  2004.  The  solar  radiation  is  superimposed  at  an 
arbitrary  scale. 


5.  Hilbert  Phase  Analysis  of  C„2,  Humidity  and 
Temperature 

A.  Phase  Locking  Between  Measurables 
To  better  understand  the  dependence  of  <I>C  upon  hu¬ 
midity  and  temperature,  we  consider  the  difference 
between  observable  phases,  as  illustrated  in  Fig.  6. 
We  define  the  following  difference  terms: 

Act  =  2  d>c  —  2  dv,  A  CH  =  2  —  2  d )#, 


Aqht  ~  2  d»c 


(2  d>T  +  2  d Ev) 
2 


(6) 


representing  the  phase  differences  between  Cn 2  and 
the  bulk  climate  parameter.  It  is  anticipated  that 
the  controlling  climate  parameter  will  be  indicated 
by  a  near-zero  difference  A.  We  reason  this  because 
if  A  is  constant  between  C2  and  a  climate  variable 
then  the  two  data  sets  must  be  phase  locked.  If  the 
mean  A  is  zero,  then  any  variations  about  zero 
should  indicate  a  synchronization  between  data  sets. 


Fig.  7.  Plots  of  Hilbert  phase  differences  between  C„2  and  solar 
radiation  for  27  and  28  March  2004.  The  solar  radiation  function  is 
superimposed  at  an  arbitrary  scale  on  each  graph.  Note  the  flatten¬ 
ing  out  of  the  phase  difference  function  during  the  daylight  hours. 


Supporting  empirical  evidence  is  found  upon  study¬ 
ing  cases  wherein  the  solar  insolation  influence  upon 
the  measured  C2  is  strong,  as  in  Fig.  7.  In  those 
plots  one  can  see  that  the  difference  between  the 
summed  Hilbert  phase  of  C  2  and  solar  radiation 
(2<I)c  _  2(hs)  flattens  out  when  the  solar  function  is 
significant. 

Motivated  by  this,  we  determined  by  linear  regres¬ 
sion  the  mean  gradients  of  all  the  A  curves,  which  we 
list  in  Table  2.  If  we  assume  that  the  only  major 
contributor  should  be  the  local  temperature  varia¬ 
tions,  Act  should  be  zero  or  near-zero  in  all  cases.  This 
turns  out  not  to  have  occurred.  In  fact,  ACT  has  no 
near-zero  value,  so  temperature  variation  cannot  be 
said  to  be  phase  locked  with  the  refractive  index 
structure  parameter. 

If  humidity  has  a  controlling  effect  in  the  visible 
and  near-infrared,  then  it  should  exhibit  a  range  of 
values  near  zero  for  ACH.  It  does  so  only  for  10  No¬ 
vember  2003  (a.m.)  and  28  March  2004  (a.m.).  It  is 
noteworthy  that  9  November  2003  shows  a  possible 
phase  lock  between  the  mean  of  humidity  and  tem¬ 
perature  with  C2.  This  seems  to  indicate  that  that 
temperature  is  vying  with  humidity  for  influence  over 
the  C2  parameter  during  that  morning. 

B.  Further  Study 

If  we  define  “dominance”  to  mean  that  the  phase 
gradient  is  the  value  closest  to  zero  of  the  differences 
under  scrutiny  then  we  find  that  humidity  is  domi- 


E  4>c  -  E  <t>H  E  <J>C  -  Z  4>t  Z  Oc  -  (  Z  ®T  +  Z  ®H  )  /  2 


Fig.  6.  Difference  plots  for  the  morning  of  10  November  2003.  The  dotted  line  is  a  linear  regression,  estimating  the  mean  phase  gradient. 
The  leftmost  graph  shows  a  phase  lock  between  C„2  and  relative  humidity. 
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Table  2.  Linear  Regression  Line  Gradients  of  the  Phase  Differences 


Date 

A, CH 

Act 

A CHT 

11/03/03  (p.m.) 

-9.7(4) 

-8.0(0) 

-8.8(7) 

11/09/03  (a.m.) 

-11.2(4) 

9.7(7) 

-0.7(4) 

11/10/03  (a.m.) 

-0.4(1) 

12.7(9) 

6.1(9) 

02/02/04  (a.m.) 

8.5(1) 

-12.3(0) 

-1.8(9) 

03/27/04  (p.m.) 

3.7(6) 

28.1(7) 

15.9(7) 

03/28/04  (a.m.) 

0.0(2) 

-3.3(1) 

-1.6(5) 

04/03/04  (a.m.) 

6.8(1) 

6.6(8) 

6.7(4) 

04/03/04  (p.m.) 

-27.3(3) 

-2.3(3) 

-14.8(3) 

nant  for  the  10  November  2003,  2  February  2004,  27 
March  2004,  and  28  March  2004  data  sets.  The  3 
November  2003  and  3  April  2004  (a.m.  and  p.m.)  data 
sets  show  that  ACT  <  ACH,  indicating  the  dominant 
effect  is  temperature,  this  being  extremely  strong  for 
the  3  April  2003  (p.m.)  data  set  and  rather  weaker  for 
the  other  two.  We  postulate  that  the  3  November  2003 
and  3  April  2004  data  sets  indicate  that  both  temper¬ 
ature  and  humidity  contribute  to  the  behavior  of  C„2, 
with  the  HPA  method  possibly  showing  the  propor¬ 
tional  contribution  to  C„2  of  each  climate  variable.  This 
is  clearly  an  area  to  be  examined  further. 

6.  Conclusion 

From  experimental  data  we  have  shown  conclusively 
that  humidity  plays  a  significant  part  in  the  visible 
and  near-infrared  measure  of  Cn2  in  a  coastal  envi¬ 
ronment.  This  is  in  qualitative  agreement  with  An¬ 
dreas’s  model  and  is  a  natural  extension  of  the  results 
from  the  mid-infrared. 

Furthermore  we  have  explored  a  new  technique, 
the  Hilbert  phase  analysis,  with  which  to  study  this 
physical  phenomenon.  In  overall  terms,  the  HPA 
method  is  in  agreement  with  the  correlogram  results. 
Phase  locking  between  data  is  an  unexpected  result 
that  needs  further  examination.  We  have  found  that 
the  HPA  technique  described  here  is  very  promising 
and  is  likely  to  provide  much  more  information  about 
the  changes  to  a  physical  system  than  traditional 
methods. 

Appendix  A.  Definitions 

The  constants  A  and  B  of  Eq.  (1)  in  the  0.36  to  3  pm 
wavelength  region  are 

A=  -10"6m1(\)(P/T2), 

B  =  4.6150  X  -  l(T6[m2(X)  -  m^X)].  (Al) 

P  and  T  are  the  ensemble  average  air  pressure  and 
temperature,  respectively.  The  functions  mx  and  m2 
are  polynomials  in  wavenumber  given  by 

6839.397  45.473 

m!(X)  =  23.7134  +  130  _  (i/\)2  +  38.9  -  (1/X)2’ 

m2(X)  =  64.8731  +  0.58058(l/\)2  -  0.0071150(l/\)4 
+  0.0008851(l/\)6,  (A2) 

with  X  units  in  micrometers. 
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